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Abstract: We investigate the radial extent of the eigenvalue distribution for Yang-Mills 
type matrix models. We show that, a three matrix Gaussian model with complex Myers 

coupling, to leading order in strong coupling is described by commuting matrices whose 

/„ \l/3 

joint eigenvalue distribution is uniform and confined to a ball of radius R = ( ^ J . We 
then study, perturbatively, a 3-component model with a pure commutator action and find 
a radial extent broadly consistent with numerical simulations. 
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1 Introduction 

Multi-matrix models arise in a wide variety of settings from matrix string theory [1], the 
IKKT model [2] (and its lower dimensional variants [3]), the BFFS and BMN models [4, 6] 
to the low energy dynamics of D-branes [7] and simple models of emergent geometry [8] 
and emergent gravity [9, 10]. More recently the BMN and BFFS models have received 
attention in the context of fast scrambling in [11]. 

In the context of numerical simulations it has been observed that the 3-component Yang 
Mills matrix model (with pure commutator action) [12] has an eigenvalue distribution that 
is parabolic with radial extent R = 2.0. It is tempting to suspect that this result should 
be derivable from first principals and this paper is an attempt in this direction. In part 
the motivation for the current paper was to attempt a theoretical estimate of the extent 
of the eigenvalue distribution in this simple model. Though we did not succeed in getting 
an exact result we succeeded in getting approximate results which are in broad agreement 
with the simulations. 

The principal results of this paper are: 

• A derivation of the uniform ball distribution in a previously studied 3-matrix model, 
some of whose observables have been solved exactly. 

• A demonstration that one can obtain the exact extent of the eigenvalue distribution 
in this 3-matrix model from an effective potential for this radial extent. 



- 1 - 



• The derivation of an effective potential to 2-loops for the pure Yang-Mills matrix 
model. 

• We derive the measure for a rotationally invariant gauge fixing which is natural if 
the matrices are approximately commuting. 

• Estimates for the radial extent in the 3-matrix Yang-Mills model, our estimate is that 
the extent of the eigenvalue distribution is R = 1.6 db 0.5. 

The structure of the paper is as follows: 

Section 2 is dedicated to the mass regulated two matrix model. This model was 
introduced in [5] and arose again in the description of the low energy dynamics of N 
L>-branes which are close together in 4-dimensional space-time [7]. The strong coupling 
properties of the model were investigated in ref. [7] and in [13] it was shown that at strong 
coupling the model is in a commuting phase with the joint eigenvalue distribution described 
by a hemisphere. 

In the first part of Section 2 we review the results of ref. [13] using a sightly different 
approach. In particular we split the matrices into their diagonal and off-diagonal elements 
and consider a general axial gauge which is equivalent to diagonalizing one of the matrices. 
We integrate out the perpendicular modes and obtain an effective action for the diagonal 
components of the matrices. Next we consider a coarse grained approximation and show 
that the longitudinal part of the diagonal modes has a parabolic eigenvalue distribution [13] . 

In the second part of Section 2 we consider an SO (3) invariant non-hermitian three 
matrix model constructed from a hermitian mass term and an anti-hermitian Myers term 
[14]. Integrating out one of the matrices one recovers the SO(2) invariant two matrix 
model considered in the first part of Section 2. However, we study the model directly in 
the three matrix model realization by splitting the matrices into their diagonal and off- 
diagonal components and consider a general axial gauge for the off-diagonal modes. Next 
we integrate out the perpendicular modes and obtain a one-loop effective action (exact in 
the axial gauge) for the longitudinal part of the diagonal components. We average over all 
possible orientations of the constant unit vector specifying the axial gauge and obtain a 
three dimensional £0(3) invariant effective action for the diagonal modes. Next we study 
the strongly coupled regime of the model and show that to leading order the ground state 
is described by a uniform joint eigenvalue distribution inside a solid ball. The hemisphere 
and parabolic distributions conjectured in ref. [13] can then be obtained by integrating 
the uniform distribution over one or two of the three coordinates respectively. Further, by 
conjecturing a uniform joint eigenvalue distribution, we obtain an effective potential for 
the radius of the distribution, which in the strong coupling limit recovers the exact radial 
extent of the distribution. It is this uniform joint eigenvalue distribution of the model 
which motivates the studies presented in Section 3 of the paper. 

Section 3 of our work explores the properties of the SO(p) invariant hermitian p- matrix 
model corresponding to a commutator squared term, i.e. a pure Yang-Mills type matrix 
model. This family of models have been studied both numerically and in a large p expansion 
in [15]. Special attention was given to the p = 3 case in ref. [16-20], including the effect 
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of finite mass and Myers terms. In particular it has been shown that, as the strength of 
the Myers coupling is varied, the ground state of the model undergoes an exotic phase 
transition from a fuzzy two sphere to a matrix phase in an exotic transition [8]. It is this 
matrix phase that we explore in Section 3. 

Section 3 contains four subsections. In the first part of the Section 3 we outline a change 
of coordinates involving splitting the matrices to diagonal and off-diagonal components. 
The change of coordinates that we consider is also equivalent to an SO(p) invariant gauge 
fixing condition for the off-diagonal elements. We calculate the corresponding Jacobian 
to second order in the off-diagonal modes. This takes into account contribution of the 
corresponding ghost modes to two loops. At tree level the ghost determinant is a generalized 
SO(p) invariant Vandermonde type determinant. 

The second part of Section 3 derives the effective action for the diagonal modes. By 
setting up a systematic perturbative expansion we obtain a two loop effective action for the 
diagonal elements of the model and under the assumption of an SO(p) invariant distribution 
for these modes within a ball of radius R we obtain an effective potential for the radial 
extent of their distribution. We show that the radial extent of these modes stabilizes in 
this approximation. For the three matrix model assuming a uniform distribution inside a 
solid ball we obtain an analytic estimate of the radius of the distribution. 

The third part of Section 3 studies the gauge dependence of our perturbative expansion. 
We consider a standard Faddeev-Popov gauge fixing procedure and compare the results 
obtained in Feynman and Landau gauges. We argue that for a generic gauge the joint 
distribution of the diagonal modes will differ from the genuine joint distribution and that 
the Landau gauge is the most suitable for studying the model. 

Section 4 of our paper is our concluding section. We summarize our results and discuss 
possible directions for future studies. 

2 Strongly coupled mass regulated two matrix model 

Let us consider the Hermitian two matrix model first considered by Hoppe [5]: 



The matrix model (2.1) arose again in the work of Kazakov et al [7] where the partition 
function for large g and a fresh derivation of Hoppe's solution for the < jj(X 2 ) > are 
given. The solution for these quantities are given in a rather complex form, which makes 
a more complete analysis of the theory (valid for large g > 1) difficult. A more detailed 
understanding of the model was initiated in ref. [13] where it was suggested that the joint 
eigenvalue distribution of the commuting saddle is given by a hemisphere distribution. 
The analysis of ref. [13] shows that at large coupling the eigenvalues of a single matrix are 
described by parabolic eigenvalue distribution. The authors then argue that this suggests 
hemisphere eigenvalue distribution for the commuting saddle of the two matrix model. 
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2.1 Hemisphere distribution 

Let us review the results of ref. [13], by following a slightly different approach which we find 
instructive. In particular we split the matrices to diagonal and off-diagonal components 
and integrate out the off-diagonal modes by imposing an axial gauge condition. To begin 
with let us define: 

Xij = Xj 5ij + a^j] Yij = Xj 5ij + o^-; Xj = (xj , x^ ); a^j = (ctjj, a^j) ■ (2-2) 
The action in (2.1) can be written as: 

N N N 

S[X,Y] = Y J X 2 i + Y, M + 9 2 l(^-^)4^ + K«")y^| 2 ; M," = l,2. (2.3) 

i=l *J=1 i , j , = 1 

Let us also consider a constant unit vector n = (n , ra 2 ) and define: 

a" = n(n.a); = (1 — nn).a . (2-4) 

Now we can use the <S?7 (iV) symmetry of the matrix model to fix the gauge: 

n.dij = . (2-5) 

Note that in this gauge the action is quadratic in the off-diagonal modes and hence the first 
loop effective action is exact! Indeed, without loss of generality we can consider n 1 ^ 
and express: 

a 1 = -%a 2 . (2.6) 
re 1 

It is then easy to check that: 

2 

a v r 1 2i n r 2 2 



or a t 



\a\a z } = T [a 2 ,a 2 ]=0. (2.7) 

n J- 



n 



This simplification is not a surprise since the gauge condition (2.5) with a constant n is 
equivalent to diagonalizing one of the matrices. The simplified action can be written as: 



N N 



S[X"] = £ [(n.x,) 2 + xf 2 ] + J2 [1 + 9 2 (n-(xi ~ ^)f] \ai 3 \ 2 + S FP . (2.8) 

i=l i,j'=l 

Here Sfp is the contribution from the Faddeev- Popov determinant resulting from the gauge 
choice (2.5). We have also split the diagonal modes to longitudinal n.x and perpendicular 
xj- components. It is easy to verify that: 



1 N 

S fp = ~2 E log(n.(fi-x,-)) 2 • (2- 



Note that if we choose a basis in which n = (1,0) or n = (0,1) the Faddeev-Popov 
determinant is just the standard Vandermonde determinant. It is now straightforward to 
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integrate out the perpendicular elements of the matrices. The resulting effective action (we 
divide by N 2 ) for the longitudinal diagonal modes is: 



1 



v 



n.Xi) 



i=l 



1 

2N 2 



N 
»,J=1 



(n.(fj - Xj)) 2 
1 + g 2 (n.(xi - xj)) 2 



(2.10) 



A few comments are in order: Note that the first loop effective action (2.10) is exact 
in this gauge. It is also valid for all g 2 . However, we know that the matrices commute 
only in the large g 2 limit. Therefore what we derived is not an exact effective action for 
the eigenvalues of the matrices but for their diagonal elements. We can define a joint 
distribution for the diagonal elements, which in the g 2 — > oo limit coincide with the joint 
eigenvalue distribution. In this spirit we consider a coarse grained description and extremize 
the following functional: 



S eS [p(x)] = I d 2 xp(x)(n.x^ 2 



d 2 xp(x) 



d 2 xd 2 x' p(x)p{x') lo£ 



(n.(x - x')) 2 



1 + g 2 (n.(x — x')) 2 



+ (2.H) 



Upon variation with respect to p we obtain 1 the equation: 

(n.(x - x')) 2 



p + (n.x) 2 = I d 2 x'p(x') lo 



1 + g 2 (n.(x — x')) 2 



(2.12) 



Now we can apply the following differential operator re.V^ to both hand sides of equation 
(2.12). The result is: 

d 2 x'p(x') 



n.x 



[ d - 

J {n.(x-x'))(l 



+ g 2 (n.{x — x')) 2 ) 



(2.13) 



This is the equation that we can use to deduce the form of the joint eigenvalue distribution 
p(x). Note that this is almost equation (5) from ref. [13]. To make the analogy complete 
we define u = n.x and choose a coordinate frame in the integral along x' in which we have 



n.x' = x . 



dx a pi(x 



Note that we have defined: 



(u-x fl )(l+g 2 (u-x f1 ) 2 ) 
Pi{x l ) = / dx 2 p(x 1 ,x 2 ) 



(2.14) 
(2.15) 



We next observe that for large g: 



- log 

2 & 



g 2 x 2 
1 + g 2 x 2 



7T 



1 



*S'(x) + 0(l/g 2 ), (2.16) 



-5{x) + 0(l/g 2 ) and - 
g x(l + g 2 x 2 ) g 



1 If we average over n here we obtain 



M + / d 3 x'p(x')4TTln(- 



1 + \ 1 + g 2 \x - x' 



- 5 - 



Substituting in equation (2.14) we obtain: 



u 



~ I dx' l p 1 (x' 1 )5\u- x' 1 ) + 0(l/g 2 ) = --p'M + 0{l/g 2 ) . (2.17) 



Next we solve equation (2.17) to leading order in g 1 and normalize f dup\{u) = 1. The 
resulting distribution is given by [13]: 

pxix 1 ) = ^(R 2 - x 12 ) , (2.18) 

with 

Now the definition of p\ from equation (2.15) and the SO (2) symmetry of the distribution 
suggests [13] the unique solution for p(x 1 ,x 2 ): 



P(^ 2 ) = ^?- (2-20) 

Equation (2.20) is the desired hemisphere distribution reported in ref. [13]. Note that the 
derivation of the two-dimensional distribution (2.20) is somewhat indirect. Indeed what 
we derived from the effective action (2.10) was the one-dimensional parabolic distribution 
(2.18). In the next subsection we will consider a three matrix model equivalent to the 
two matrix model (2.1) and prove directly that at strong coupling the corresponding three 
dimensional eigenvalue distribution is an uniform distribution inside a solid ball. The two 
dimensional hemisphere eigenvalue distribution can be obtained by integrating out one of 
the eigenvalues. 

2.2 Three matrix model realization and uniform distribution 

Let us now consider the model (also originally introduced parenthetically by Hoppe [5] 
page 73 and further discussed in [7] and [13]: 



j VX VYVZe- N ^ x2+Y2+z2 - ia ^ Y W . (2.21) 



It is easy to verify that if one integrates out the Z matrix and defines g 2 = (ia) 2 /4 
one recovers the two matrix model (2.1). This suggests that the model (2.21) should be 
as solvable as the two matrix model. Note also that there is a global SO (3) symmetry 
rotating the X, Y and Z matrices. We find it instructive to analyze the model in the spirit 
described in section 2.1. To begin with let us define: 

Xij = x^5ij-\-(i^j', Yij = x^5ij-\-a,^j', Zij = x^5ij-\-ci^j', Xi = (xj,Xj,x ); Ojj = {ci^j , cl^j , Q>ij) ■ 

(2.22) 

Next we consider a constant unit vector n = (n 1 ,^ 2 ,?! 3 ), define: 



and impose the axial gauge n.a = 0. The action in (2.21) can then be written as: 

N N 



Six, a] 



£[(«■ 

i = l 



*ij'=i 



. Ot 



-i-e^pnPin.Aij) 



(2.24) 



where we have defined Aij = xi — xj. Note that there is no term cubic in a in the action 
(2.24), because of the axial gauge. This suggests that the first loop effective action is exact. 
Now we proceed as in section 2.1 and integrate the perpendicular matrix elements. One 
can show that the resulting effective action is given by: 



1 N 1 N 

S eS [(n.x)] = ^^(n-Xi) 2 -^ Yl lo ^ 



i=l 



g 2 (n.(xj - xj)) 2 
1 +g 2 (n.(xi - Xj))' 



(N-l) n , 
+ -^T^9 2 , (2-25) 



where g 2 = (ia) 2 /4. Next we consider a coarse grained approximation an vary the corre- 
sponding distribution function p to obtain the equation: 



p + (n.x Y = I d 3 x'p(^') lo 



g 2 (n.(x — x')) 2 
1 + g 2 (n.{x — x')) 2 



(2.26) 



Note that equation (2.26) is valid for any choice of n. Next we average over all possible 
directions that n can take with a uniform weight. Or equivalently average over the unit 
two-sphere. It is easy to show that: 



— / dtt 2 {n.x) 2 = \x 2 
4vr J d 



(2.27) 



The right-hand side of equation (2.26) requires a bit more careful analysis. One can show 
that: 



J(g\x-x'\) := i- / 



g 2 (n.(x - x')) 2 
1 + g 2 (n.(x — x')) 2 



g\x-x'\ 



2g\x — x'\ 



di] log 



rf 



-g\x—x' 



-2 arctan(g(|x — x'\) 
g\x — x'\ 



+ log 



g 2 {x — x') 2 



1 + g 2 (x — x') 2 



7T 1 



g \x — x' 



To leading order in 1/g equation (2.26) becomes: 



1 ^2 
p + -x 



-- f d 3 x'p(x')—±— 

gj \x-x'\ 



l + r] 2 
(2.28) 

(2.29) 



Next we apply the operator of Laplacian A x on both hand-sides of equation (2.29) to 
obtain: 



4vr^ 
9 



Therefore we conclude that: 



d 6 x'p{x')5{x-x*) 



p(x) = — j = const 



4vr^ 
9 



p(x) 



(2.30) 
(2.31) 
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and hence the ground state at strong coupling corresponds to an uniform eigenvalue dis- 
tribution inside a ball. In order to estimate the radius of the eigenvalue distribution R we 
use the normalization of p(x): 

J d 3 x P {x) = y r3 P = 1 • ( 2 - 32 ) 
From equations (2.31) and (2.32) we obtain: 

«=(!) 1/3 , (2.33) 

which is exactly the radius obtained in ref. [13] reported in equation (2.19). It is also 
straightforward to obtain the hemisphere distribution (2.20). Indeed: 



vV-xi 2 -^ 2 _ _ 

P2(X\X*) = J dx*p{x\x\x*) = 3 — • ( 2 - 34 ) 

-y/R 2 -x l2 -x 22 

Some additional comments are worthwhile here. First one can push the analysis further 
by observing that if we don't assume large coupling in stead of (2.29) we obtain 

p + ^x 2 = J d 3 x'p(x')J(g\x - x'\). (2.35) 

which upon acting with the Laplacian and noting 

V 2 J(x) = — (2.36) 

(1 + g 2 \x — x'| 2 )|x — x'\ 2 

we obtain the integral equation 

1 = / dV / ( ? - _ (2.37) 

J (1 + g 2 \x — x'| 2 )|x — x'\ 2 

whose large g behaviour is 

1 - 7T5 ^~ X/l) +0(l/g 2 ) = 2 -^5(x-x') + 0(l/g 2 ) (2.38) 



/|2 



(1 + g 2 \x — x'\ 2 )\x — x'\ 2 %g \x — x 

So to leading order we obtain p{x) = ■^ 7 6(R — r). 

Given that the eigenvalue distribution is concentrated in the interior of a ball of radius 
R we can further deduce that this radius is determined by 

9 2 = ^{ 9 R? + --- and v= 9 ^ f p(x)x*= { ^f + ... (2.39) 
we find that to leading order in large g the observable 
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in agreement with [7]. 

Further: assuming a constant eigenvalue distribution within a sphere of radius R i.e. 
P( r ) = 4§jip and usin S 

<Px<Px' P £ )p W = -4 (2.41) 
\x — x'\ 5R 

we obtain 

^ (fi) = T + |k (2 - 42) 

Varying with respect to R we find that R? = || in agreement with the exact expression 
(2.19). In the next section we attempt to use such an effective potential (derived perturba- 
tively) to estimate the extent of the eigenvalue distribution in a three matrix model with 
pure commutator action. 

3 The p > 3 matrix model 

In this section we consider the p-matrix model: 

S[X}=NTr(-±[X»,X»A; l x,p = l...p; (3.1) 

where X a are hermitian N x N matrices. The partition function is given by: 

Z = J dX a e~ s ^ = J VXe N \ T ^^\ (3.2) 

Note that in addition to the SU(N) gauge invariance the model (3.2) has a global SO(p) 
symmetry transforming the matrices X^ . We are interested in the eigenvalue distribution 
of one of the matrices and in particular in the extent R of the eigenvalue distribution in 
the large N limit. We find it convenient to split the degrees of freedom to diagonal and 
off-diagonal contributions: 

= + (3-3) 
In terms of the new variables the action in (3.2) can be written as: 

S[X]=ivi£ 14/^11^ , (3.4) 

where we have defined: 

nJ = Ag/|^|; Hf/ = ^ - n^; . (3.5) 

A standard way to proceed would be to integrate out the off-diagonal degrees of freedom 
afj and obtain an effective action for the diagonal components xf. Note that the quadratic 
term in in equation (3.4) is proportional to a projector and hence cannot be directly 
inverted. A gauge fixing is required. We find it natural to work in a gauge in which 
the longitudinal modes are removed, more precisely we impose the gauge fixing condition 
fiij.dij = 0. In the next subsection we briefly present the change of variables necessary 
to implement our gauge fixing condition. We refer the reader to Appendix A for a more 
detailed calculation of the corresponding Jacobian. 
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3.1 Gauge fixing 

Our goal is to perform a change of coordinates which is manifestly SO(p) invariant and con- 
venient in calculating quantum corrections to the effective potential governing the ground 
state of the theory. 

Let us consider a set of p— 1 hermitian matrices aj- 1 ' (m = 1 . . . p— 1; i,j = l... N) with 
vanishing diagonal components (a™ -1 " = 0). A slightly more general way to parametrize 
this set of matrices is to consider a set of p linearly dependent matrices a^- (fi = 1 . . . p) 
satisfying: 

p 

^V.a^ = n.a L = , (3.6) 
n=i 

where n is a constant p— dimensional unit vector and we have suppressed the indices 
Next let us consider any set of TV 2 — TV orthogonal matrices Rij £ SO(p) satisfying: 

Rij.n = fiij for i>j; (3-7) 
Rij = Rji for i < j; ■ 

Clearly such a set of matrices always exists. Now we define the hermitian matrices: 

Sij = Rij-ajj for i j and an = 0; . (3-8) 

The matrices linearly dependent and satisfy the properties: 

fiij.dij = (3.9) 
aij.aj i = aj- j .aj i (3.10) 

Next we define the change of coordinates: 

X 11 = U{x l * + a l *)U- 1 ; U G SU(N); . (3.11) 

Note that on the left hand-side of equation (3.11) we have a set of p N x N hermitian 
matrices spanning a pN 2 dimensional linear space. On the other side by construction there 
are only p — 1 linearly independent matrices labeled by and hence the dimension of the 
linear space spanned by + a M is equal to (p — l)(iV 2 — N) + pN. This suggests that 
in order to have a well defined change of coordinates in equation (3.11) we need iV 2 — iV 
degrees of freedom, parameterizing the orbit of + under the SU(N) group. However 
a general element of SU (N) has iV 2 — 1 degrees of freedom. The N — 1 degrees of freedom 
that are left out correspond to the stabilizer of the action of SU (N) and as one can check 
are generated by a Cartan subalgebra of su(N). 

It is a straightforward exercise to compute the corresponding Jacobian. The expression 
that one obtains is given by (we refer the reader to Appendix A for a detailed computation) : 



J 




d i°j +Y ij 



det 



56 r 



(3.12) 
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where Y-J 11 is given by: 

Y\T = Tr-A***? - + Uij ^ j f 2 ml .{5? -Si- 5? + 4) . (3.13) 

I ^ij I I ^ij I 

Note that the last determinant in equation (3.12) is the Haar measure of SU(N). Now we 
can write down the measure in the path integral (3.2) in terms of the new variables. To 
second order in the off-diagonal elements Sij we have the expression: 



1 - 2 3ij.Dij.Sji + 0(a 3> 



(3.14) 

where DC/ is the Haar measure of SU(N) and Dij is given by: 

nM v _ n f/ 1 \ - n ii n ji + n »7 n jZ /„ -, ,s 

« + ^ EPS ' ( } 

Note that without loss of generality one can take n = e p . Note also that By = Rij.Sjj. 
Our next goal is to develop a systematic perturbative procedure to integrate out the off- 
diagonal degrees of freedom ajj . 

3.2 The effective action to two loops and stabilization 

In this subsection we will develop a perturbative technique to integrate out the off diagonal 
modes Sjj and calculate the semi-classical correction to the effective action for the diagonal 
modes Xj. 

Our experience with the two matrix model from section 2 and in particular its three 
matrix model realization suggests that the ground state of the theory corresponds to a 
uniform joint eigenvalue distribution inside a ball of radius R. Furthermore numerical 
studies of the three matrix model [12] give radius R ~ 2.0 > 1. This suggests a perturbative 
expansion in powers of 1/R may prove useful. In order to proceed systematically we first 
rescale our variables Xj and Sij in the following way: 

x = Rxf, 3 = RSij] . (3.16) 

Next we write the action (3.4) as: 

S[k, = NR 4 £ (-\\k l3 \%.~3^ + A^a", a"]* + \[a", ti%[d», a v ] jt ) (3.17) 

Note that from field theory point of view the parameter l/NR 4 can be interpreted as a 
loop counting parameter. Next we define the correlation function: 



(O) = * — '- -± 3 -L = -j— j- . (3.18) 

SUidxi (IL* IU>j (dP-^dp-^aA e-N&^h^l 2 ^ 
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For the propagator of we obtain: 

NR 4 



\ a ij a lm/0 - 7754 ~ <>i Oj . (d.iyj 

l A U'| 2 



In deriving (3.19) one could use a frame in which n = e p . Note that the cubic and quartic 
contributions to the action (3.17) as well as the Jacobian in the measure (5.13) depend 
on through the relation a^j = Rij.a^j. This is why it is convenient to calculate the 
two-point function: 

«f a im)o = ^4 ~} ^1?^) , (3-20) 
I Ajj I 

where we have used that: 

Rtf.(5^' - n»'n v ').R\f = ST - n%n% = 11% . (3.21) 

We now have all the machinery required for a perturbative calculation in powers of l/NR 4 . 
To leading order we obtain the following first loop effective action for the diagonal modes: 

SS ( R i = l(P ~ 2 ) n2 ~ 2(P " l)N + p] log R + Yl lo § l^^'l 2 ' ( 3 - 22 ) 

As one can see from the first term in equation (3.22) at large N and for p > 2 at one loop 
the effective action gives an attractive potential V c s(R) and is therefore not sufficient to 
stabilize the radius of the distribution. On the other hand the second loop corrections has 
an overall factor of 1/NR and could balance the log-R attractive potential in (3.22). 

At second loop the effective action has contributions from the cubic and quartic vertices 
in (3.17) as well as from the quadratic term in the measure (5.13) ("ghost's contribution"). 
The corresponding correlation functions are: ((iVi^A^-aJJa^, a u ]ji) 2 )o, {N R 4 \tr[a^ , a?] 2 )o 
and {—2Y^i jO-ij.Dij.(iji)o. Using Wick contractions and the two-point function (3.20) we 
can calculate the second loop contribution. After somewhat tedious but straightforward 
calculations we obtain: 

/f/VR4AM~™ ~ v] ^ _ 1 y i( 4 P ~ 6) sin 2 djji + sin 2 hij } 

< W «« W an\ = ^j: VP-W : *l + *°'««} , (3 . 24 ) 
4 2NR^ f-ri A 2 A 2 

sr~ v — 1 v^. 2 sin 2 9, ■ ~i „ / 1 \ 

(-2yj 3i3 ,A 3 , 3ji >0 = ^ £ -gZ* +0 [jfj ■ (3.25) 

M hJ^l ^il^jl 

Were the angles 6i t ij,9iji are defined via cos O^ji = riij.nn and the last term in equation 
(3.25) corresponds to non-planar diagrams contribution subleading in the large ./V limit. 
For the total second loop contribution to the effective potential we obtain: 

s$ (R, a) = {J ^4 V p- 1 : 4 !^ 9 ^ . (3.26) 

1,3 & ^il^jl 
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The full large N effective action for to one loop is then 

Since the ground state of the theory is in a commuting phase we can trade the discrete 
sums in equations (3. 22), (3. 26) for integrals over a joint eigenvalue distribution p(x) via: 



where the integral is over a ball B p of unit radius. In the large N limit we obtain the 
following second loop effective potential for the radius of the joint eigenvalue distribution 
R: 

V eS (R) = (p - 2)N 2 (logi? + , (3.29) 

where we have defined: 

#(P) = \l I I &xd^dVzp(x)p(y)p(z) V - 1 ~**f |f (3.30) 
^JbpJbpJbp {x — z) z (y — z) z 

and x ,yz is the analog of diji defined via: 

(x-y).(x- z) 

costi xyz = — = ~ = — . (3.31) 

\(x - y)\\(x - z)\ 

We can now estimate the radius of the joint eigenvalue distribution corresponding to the 
ground state of the theory. Indeed minimizing (3.29) we obtain: 

R p = (4#(p)) 1 /4 ; . (3.32) 

Note that for p > 5 we have > (the integrant in equation (3.30) is non-negative) 

and the ground state stabilizes at the finite radius estimated in equation (3.32). For p = 3 
we evaluated analytically the integral in equation (3.30), assuming uniform joint eigenvalue 
distribution p(x) = const (look at Appendix B for more details). The resulting radius is: 



1/4 

R 3 = I 9 - ^ ) » 1.323; . (3.33) 



3 .2 



Equations (3.29)-(3.32) contain the main result of our perturbative calculation. In the next 
section we discuss the gauge dependence of our expression for the radius of the distribution 
R p in equation (3.32). 

3.3 Gauge dependence 

In subsection 3.1 we outlined a change of coordinates that was equivalent to introducing the 
gauge fiij.dij = 0. Alternatively we could have used a standard Faddeev- Popov techniques 
to fix our gauge. Let us consider the gauge condition: 

fij = Aij.Sij = 0; . (3.34) 
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The corresponding Faddeev-Popov determinant is given by: 

llm i 

\i\ij\-aex,\\o\oy + r 

i>j 

where: 



A FP = 11 |A^| 2 det||^f + Y'$\\ , (3.35) 



I ^ij I I ^ij I 

Notice that the gauge condition fV. = Aij.Sij — rjij = would result to the same Faddeev- 
Popov determinant (3.36). Now integrating over the family of gauge functions with 
weight exp(— |r/jj| 2 /2^) would modify the action (3.4) to: 

S[X,t] = A^IA./a^Iir + " Nj^^^ary - N±tr[a»,a»f , 

(3.37) 

Next we can go through the steps considered in section 3.2, namely rescale with the radius of 
the joint eigenvalue distribution R as in equation (3.16) and set up perturbative calculation 
in powers of 1/R. One can show that the first loop effective action is still given by 
equation (3.22) and is thus gauge independent. However the two-point function (3.20) (the 
propagator for a^) is modified to: 

Wtfmk = pi^ff + HK)C4 • (3-38) 

| Ay | 

Note that the result from equation (3.20) corresponds to the choice £ = (Landau gauge). 
Let us consider the choice £ = 1 (Feynman gauge) and calculate the second loop contribu- 
tion to the effective action. Going through the same steps as in section 3.2 we obtain the 
analog of equation (3.29): 

V cS (R) = (p-2)N 2 ^ogR+#j^j , (3.39) 

where is given by: 

#(P) = V -^r\ I I dVy<F~zp&p($)pCz) 1 (3.40) 

z JbpJbpJbp {x — z) z (y — z) z 

This results to the radius: 

R' p = (4#(p)) 1 /4 ; . (3.41) 

It is clear from equation (3.40) that is positive for p > 3. For p = 3 we have evaluated 
analytically #(3) (look at Appendix B). The corresponding radius is: 

R 3=[l + f^J « 1-797 . (3.42) 

Apparently the results obtained in Landau and Feynman gauges differ. In order to 
address the issue of gauge dependence let us focus on a particular representative of the 
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family of gauge conditions fV. = Aij.dij — rjij = 0. Note that the change of coordinates 
(3.11) considered in section 3.2 implements the rjij = case. One can show that the gauge 
condition for general rjij can be implemented along the lines of section 3.2 via the following 
modified change of coordinates: 



n 



X» = UWxTSy + RgaZf + -^-TjijWU- 1 ; U e SU(N); . (3.43) 



A 



Let us suppose that the theory have settled in its ground state which is a commuting phase. 
There should exist unitary matrix V £ SU(N) which simultaneously diagonalizes the X^ 
matrices and hence we can write: 

u 

n 

V~ l \»V = UWxfSij + R^a^ + -f-MjWU" 1 , (3.44) 

where A M is a diagonal matrix. Now if we square equation (3.44), take a trace over the 
gauge indices and sum over \i we obtain: 

E 3 = E *1 + E 1 4 1 2 + E M ^ E -1 + E S ■ ( 3 • 45 > 

i % ij ij I *J I j ij \ % 3\ 

Next we define average radii of the distribution and r x via: 



-a = ^E^ ^ = ^E-1; (^6) 



and learn that: 



Therefore the average radius of the eigenvalue distribution r\ always differs from the aver- 
age radius of the distribution of the diagonal modes r x , unless rjij = or the eigenvalues are 
infinitely spread in which case there is no well defined average radius. This could explain 
why the gauge fixing procedure outlined above, which involved averaging over all possible 
values of rjij failed to produce a gauge independent answer for the radius of the eigenvalue 
distribution. These consideration suggests that the gauge rjij = should be optimal for 
describing the almost commuting theory. 

Alternatively one could take the point of view that both gauge choices are equally valid 
and describe different approximations to the true result, they only differ due to the intrinsic 
errors in a perturbative calculation. If we take this point of view we can use the difference 
to estimate the errors in our estimate of R. If we do this we conclude that R ~ 1.6 ± 0.5 
which is in reasonable agreement with the numerical results. 

4 Discussion 

In this paper we have followed two treads, in the first we investigated the 3-matrix model of 
[5, 7] and find that in the large g limit the 3-matrices commute and have a joint eigenvalue 
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/„ \l/3 

distribution given by the uniform distribution within a ball of radius R = ( ^ J . We 
show that a simple effective potential for the radius of the distribution reproduces the exact 
result. 

Encouraged by the success of this effective potential calculation we develop an effective 
potential for the radius of the p-component Yang-Mills matrix model to two loops. We have 
done this by deriving an effective potential for the diagonal modes while preserving SO(3) 
invariance and then assuming that these modes are uniformly distributed. The direct 
analog of the computations in earlier sections would be a two loop computation in the 
axial gauge (where one of the matrices is diagonalized) . Unfortunately this gauge choice 
leads to infrared divergences at two loops and so we have not pursued this option. 

We found that it is necessary to go to two loops as at one loop the effective potential 
is not stable since there is no classical potential and the one loop term gives an attractive 
potential which is a rotationally invariant version of the Vandermonde determinant. The 
eigenvalue repulsion arises at two loops and gives a hard core potential. It is easy to 
see that higher order terms give inverse higher powers of R and our two loop estimate can 
only be a very rough approximation. 

Our estimate for the radius is unfortunately gauge dependent with R = 1.323 in the 
Landau gauge and R = 1.797 in the Feynman gauge. It is reasonable to assume that the 
difference between these is an indication of the errors in the method which would indicate 
that perhaps a reasonable estimate can be obtained by averaging the two and taking the 
difference as an indication of the error yielding the prediction R = 1.6±0.5. An alternative 
approach pursued by Hotta, Nishimura and Tsuchiya [15] examined similar questions in 
the general Yang-Mills p-matrix model 2 analysis of some observables. If we take their result 

<^)>=/f(l + -U...) (4.1) 

and assume that this is valid for p = 3 together with the assumption that the eigenvalue 
distribution of a single matrix is parabolic of extent R (which is consistent with a uniform 
joint distribution for commuting matrices) then < |^(X 2 ) >= 3-^- and using (4.1) we 
obtain the estimate R = 1.68, which is surprisingly close to the estimate we obtain above. 
In this case one can also attempt an estimate of the error by noting that if instead of (4.1) 

we use < %(X%) >= a/Itttxt + • • • ) which has the 

same leading large p expansion we 

obtain R = 1.83 and therefore estimate the error (within the assumption of a parabolic 
distribution) that R = 1.75 ±.15 which is in surprisingly good agreement with Monte Carlo 
simulations [12] which give a value for R = 2.0 

We conclude that, though the random matrices are not commuting, a useful approxi- 
mation is to take the background formed by these fluctuating matrices as that commuting 
matrices whose joint eigenvalue distribution is approximately uniform within a ball of ra- 



2 Hotta et al [15] also developed calculated a two loop effective action for eigenvalues of the matrices X M 
by integrating out the U (N) transformations that diagonalize the matrices. This yields a non-rotationally 
invariant effective action, which has does not lead to a stable effective potential for the radial extent of the 
eigenvalues. 
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dius R. This gives background reasonable agreement with the numerical work and serves 
as a reasonable starting point for further work. 



5 Appendix 
5.1 Appendix A 

In this subsection we outline the calculation of the Jacobian (3.12). Let us begin by 
differentiating equation (3.11). We obtain: 

(U^dXUjij = dxidij + dRij.ajj + Rij.dajj - \ Aij\ftij0ij - [a,6]ij , (5.1) 

where we have defined the Maurer-Cartan form 9 = U~ 1 dU. Next we define the tetrads E 
via: 

Ei = (U^dXUU Eij = (U^dXU)^ for i j= j; . (5.2) 
In matrix form we have: 



\E\\ = 





dx k 




duim 




16} 





S[a,6}u 
Sui m 


4 


IT dRij a • 
iil .r dx k 


1616] 


l 3' Suim 


4 


dfiij -> 
dx k -"y 





IA 1 I n 



where IT* = 1 — HijUij and we have split the off-diagonal tetrads Eij into: 



/••',, 1 1 ;,./••„: and £ 



We have also used that: 



(5.3) 

dx k ' dx k '" LJ " ' 

Now the Jacobian that we are interested J is given by the determinant of \\E\\. It is an 
easy exercise to show that the Jacobian is given by: 



J = det 



A-- 1 3 



lJ $ui m %r 5u lm 



5[a, 9]ij dfjjj _, 5[a, 9} kk 

^ dx k ' 13 5ui m 



(5.5) 



-17- 



One can show that the determinant in equation (5.5) factorizes: 



J = det 



<9r% _ 5[a^,6] kk 

• On' 



det 



59 r 



5u ln 



(5.6) 



The last determinant in equation (5.6) produces the correct measure of SU(N) it is the 
first part J' that we are really interested in. One can verify that: 



J'= (ni4il]det 



d s;m i \rlra 
°i 6 j + Y ij 



(5.7) 



where is given by: 



<^jl - II jo .fl m / 



- 'J 



IA 



IA--I 2 



where we have used that: 



56i rn 
Now using that: 



- (aizc)j - a mj dj), —7 — t^— (Oj - dj); ajj.llj-,- — ajj, 

|Ajj| 



(5.8) 



(5.9) 



det 



d i d j +Y ij 



exp <| try - -try 2 + 0(Y 3 ) 



One can easily verify that: 



(5.10) 



J' = exp < 



^log|Ajj| - 2^J%.Aj-Qji + Q(e 



(5.11) 



where -Djj is given by: 



2 + 4 

Therefore our final expression for the measure is: 



|Ajfc||Ajfc| 



(5.12) 



DU Y[ dxi Yl (dP^aijdP^aji) exp I log |Ajy| - 2 ^ ^.Dij-aji + 0(c 



i>j 



i -J 



'■J 



(5.13) 



5.2 Appendix B 

In this subsection we will provide details about the analytic evaluation of the quantities 
#(p) and #(p) defined in equations (3.30) and (3.40) respectively. Note that one can write: 



#(P) = + 2C P ; #(p) = V -^A P - , 



(5.14) 
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where: 

A p = [ I f tFxePy<?zp&ptf)p@) ~ * - (5.15) 
JbpJbpJbp (x — z) z {y — z) z 

C p = I [ [ d p xd^zp(x)p(y)p(z) °°f ^ (5.16) 
JbpJbpJbp (x — z) z \y — zy 

Let us consider first the quantity A p . Note that the integrals along x and y in equation 
(5.15) factorize and one can write: 

2^/ 2 



BP 



r(p/2) 



A p = / dPzp(f)Q p (|i|) 2 = I dzz^ x p{z)Q p {zf , (5.17) 

where: 



Q p (z) = [ dPxp{x)-J^- (5.18) 

J BP \X — Z) z 

and we have used that the distribution is SO{p) symmetric (namely p(x) = p(x)). Note 
that for a uniform distribution we have p(x) = pT(p/2) /2nP/ 2 and we have: 

i 

A p =p j dzz p - l Q p {zf . (5.19) 



One can show that for p = 3 and a uniform distribution Qs(z) is given by: 

3 



Qs(z) 



1 - z 2 
1 H Arctanh(i) 

z 



(5.20) 



and 

l 

A 3 = 3 y ^5 2 Q 3 (5) 2 = ^ (15 + 2vr 2 ) . (5.21) 
o 

Let us now focus on the quantity C p . After using equation (3.31) and the SO(p) symmetry 
of the eigenvalue distribution, one can show that: 

Cp = -A p -I [ [ d?xd*yd?zp{x')pQ)p& £ (5-22) 
4 JbpJbpJbp (x — z) z (y — zp 

3 27r p / 2 f 

= + y ^F-V^Q;^)^^) , (5.23) 



o 

where 



%(z) = -l d p xp(x)log\(x--~z)\ . (5.24) 

* J Bp 

For an uniform distribution p(z) we have: 



i 



C p = ^A p + pJ dzzP^Q'pW&piz) . (5.25) 
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One can show that for p = 3 and a uniform distribution $3(5) is given by: 

5 log 1 1 + 

(5.26) 



$3(5) = — [_34z + 65 3 + 3{z - l) 3 (3 + 5)log|l -51 - 3(5 - 3)(1 + 5) 3 log 1 1 +511 
96z 



and one can calculate C3: 

c ^l( i5+ H- <5 - 27) 

Now one can substitute the results from equations (5.21) and (5.27) into equation (5.14) 
to obtain: 

#(3) = 1(15 - vr 2 ); #(3) = 1(15 + 2vr 2 ); , (5.28) 
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